Pe 


arX 


0707.0208v1 [astro-ph] 2 Jul 2007 


1V 


Foam-like structure of the Universe 
A.A. Kirillov and D. Turaev 


Abstract 


On the quantum stage spacetime had the foam-like structure. When 
the Universe cools, the foam structure tempers and does not disappear. 
We show that effects caused by the foamed structure mimic very well 
the observed Dark Matter phenomena. Moreover, we show that in a 
foamed space photons undergo a chaotic scattering and together with 
every discrete source of radiation we should observe a diffuse halo. We 
show that the distribution of the diffuse halo of radiation around a point- 
like source repeats exactly the distribution of dark matter around the 
same source, i.e. the DM halos are sources of the diffuse radiation. 


1 Introduction 


Wheeler pointed out that at Planck scales topology of spacetime undergoes 
quantum fluctuations [I5]. In the present Universe such fluctuations carry vir- 
tual character and do not lead to detectable topology changes. However, in 
the past, the Universe went through the quantum stage when the temperature 
exceeded the Planckian value and the fluctuations were strong enough to form a 
non-trivial topological structure of space. In other words, the very early, quan- 
tum stage Universe had to have a foam-like structure. During the cosmological 
expansion, the Universe cools, quantum gravity processes stop, and the topo- 
logical structure of space freezes. There is no obvious reason why the resulting 
topology has to be exactly that of R? — relics of the quantum stage foam might 
very well survive, thus creating a certain distribution of wormholes in space. In 
the present paper we show that the whole variety of the observed Dark Mat- 
ter (DM) phenomena admits a straightforward interpretation in terms of the 
foam-like topological structure of space. Moreover, the specific properties of the 
foam that are read from the observed DM distribution coincide with those that 
are derived theoretically from very basic physical principles: we show that the 
actual distribution of DM sources corresponds to the ground state of the linear 
field theory on the foamed space. 

An arbitrary non-trivial topology of space can be described as follows. Given 
a Riemanian 3D manifold M, we take a point O in it and issue geodesics from O 
in every direction. Then points in M can be labeled by the distance from O and 
by the direction of the corresponding geodesic. In other words, for an observer 
at O the space will look as R? (endowed with a metric lifted from M). Given 
a point P € M, there may exist many homotopically non-equivalent geodesics 
connecting O and P. Thus, the point P will have many images in R. The 
observer might determine the topology of M by noticing that in the observed 
space R? there is a fundamental domain D such that every radiation or gravity 
source in D has a number of copies outside D. The actual manifold M is then 
obtained by identifying the copies. In this way, we may describe the topology 
of space M by indicating for each point r € R? the set of its copies E(r), i.e. 


the set of points that are images of the same point in M. Most of the time, we 
will simply speak about the images of points in R3, without referring to M. 

Note that an observer ignorant of the actual topological structure of M will 
greatly overestimate the density of matter (as all the gravity sources outside 
the fundamental domain D are fictitious — each of them is just an image of 
some point in D seen from another direction). However, one cannot immedi- 
ately apply the above picture to the explanation of DM effects: the Dark Matter 
emerges on galaxy scales while we do not see multiple images of galaxies densely 
filling the sky. Our idea that allows to link the observed DM effects with the 
topological structure of space is that the fundamental domain may be of such 
distorted shape that the direct recovery of the actual topology of space by de- 
tecting images of sources could be impossible. Indeed, the non-trivial topology 
at present is a remnant of quantum fluctuations at the very early Universe, and 
the randomness built in the structure of the original quantum foam can survive 
the cosmological expansion. Namely, at the quantum stage the state of the 
Universe was described by a wave function defined on the space of Riemanian 
3D manifolds. Once quantum gravity processes stop, the further evolution of 
the wave function was governed by the cosmological expansion only. It is highly 
unlikely that the expansion could led to a complete reduction of the wave func- 
tion, i.e. to singling out one definite topological structure of the Universe. In 
other words, if at the end of quantum gravity era the Universe was not in a 
pure quantum state, it is not in such a state now. One cannot, therefore, speak 
about a definite topological structure of space, i.e. assign a definite set E(r) 
of images to every point r € R3. A point r’ € R can be an image of r with 
a certain probability only, hence instead of a discrete set of images, a smooth 
halo of images of every single point appears. 

Even if we want to believe that a definite (classical) topological structure has 
happened to emerge out of the quantum foam, the randomness of this structure 
will persist: the wormholes which remained as the quantum foam tempered 
will be randomly cast in space. Moreover, we recall that a typical wormhole 
is obtained as follows: the interior of two remote spheres is removed from R 
and then the surfaces of the spheres are glued together}. Such wormhole works 
like a conjugated couple of convex (spherical) mirrors, therefore a parallel beam 
of geodesics diverges after passing through the wormhole. Thus, if we place 
spherical wormholes randomly in R?, the flow of geodesics that pass through a 
large number of the wormholes will have a mixing property (like the flow of Sinai 
billiard, or of Lorenz gas). For a point-like source for radiation or gravity, it 
means that some portion of photons/gravitons will be scattered by the spherical 
wormbholes, which will create a specific smooth halo around every single source. 

In any case, no matter what is the exact origin of the randomness of the 
topological structure of space, one can take such random structure into account 
by introducing a certain measure on the space of all Riemanian 3D-manifolds 
M. The observed topological or metric properties of space are then obtained 


tone can imagine a more general construction as well, where a pair of more complicated 


two-dimensional surfaces replaces the spheres 


by averaging over this measure. Thus, for example, an individual manifold M 
is defined by specifying, for any point r’ € R? the set E(r’) of its images (the 
points in R? that represent the same point of M). Averaging over all manifolds 
M, gives a distribution K(r,r’) of the images of r’: 


K(r,r’) = 6(r =r’) + O(r, 1’), (1) 


where the first term corresponds to the point r’ itself, while 6(r,r’) is a certain 
smooth distribution of additional images of r’; namely, in the neighborhood of 
a point r of volume d?r there is (on average) b(r,r’)d?r images of r’. 

It means that a single particle of matter at the point r’ is always accompanied 
by a smooth density b(r,r’) of exactly the same matter. This halo does not 
necessarily emit enough light to be identified, but it will always contribute to 
gravity. Thus, if the halo is not seen, it is detected by an anomalous behavior 
of the gravitation potential of the point-source. Such anomalous behavior is 
indeed universally observed starting with the galaxy scales, and constitutes the 
DM phenomenon. The existence of a quite rigid dependence between the density 
of luminous matter (LM) and the density of DM is a well-known observational 
fact. This fact allows us to interpret the DM phenomenon as an indication of 
the random topological structure of space, with formula giving 


ppm(r) = [oe r)prm(r) dr. (2) 
In fact, the simple law 
blr, r) ~w |r—r'|-? at |r -r'| > Ro (3) 


(where Ro is the galaxy scale) provides quite accurate description of all known 
DM effects. In particular, it allows to recover the whole variety of observed 
galaxy rotation curves [6]. It is also consistent with the observed fractal struc- 
ture of the distribution of matter on large scales [4 [3] [0]. 

Note that relations (2),(8) give a good description for the observed DM phe- 
nomena, independently of a theoretical interpretation [4] [6]. We will, however, 
show that in our picture where b(r,r’) is an averaged characteristic of the topo- 
logical structure of space, empirical law acquires a basic physical meaning. 

It is also important that in our interpretation the DM halo is not actually 
dark. The image r of a point r’ represents the same physical point, just seen 
from another direction. Therefore, if the source of gravity at r’ is also a source 
of radiation, all its images in the halo will be luminous too. However, the 
halo radiation has a diffuse character and the brightness is very low (the halo 
radiates a reflected light, in a sense). In observations, relating the halo radiation 
to a particular point source could be a very difficult task. In fact, the presence 
of a significant diffuse component in cosmic radiation, unidentified with any 
particular source, is well known [12]. Usually, the observed diffuse halos in 
galaxies are attributed to reflection from dust, and the general diffuse component 
is assumed to originate from very fade and remote galaxies, but it has never 


been related to DM halos. However, it was very convincingly demonstrated in 
that the observed DM/LM ratio within the intracluster gas clouds is much 
less than that for galaxies. This observation gives a strong argument in support 
of our theory of DM effects: while for small and bright sources (galaxies) the 
luminosity of the halo is filtered out by the observer and the halo appears to 
be dark, for the extended radiation sources (cluster size plasma clouds) the 
diffuse halo radiation comes from the same region of space and is automatically 
accounted in the total luminosity of the cloud. 

Indeed, we show below that the intensity of sources of radiation renormalizes 
according to the following law: 


Liotal (r) = Lsource (r) + Ihalo (r), (4) 


where 
Thato(r) = fe, it Looe? Or, (5) 


with the same b(r,r’) as in (2). Therefore, in our picture, the luminosity of the 
DM is always proportional to its density. The gravitating halos of discrete light 
sources in the sky only appear to be dark, because of their diffuse character. 

From the physical standpoint the foamed space is a porous system. It means 
that the coordinate volume, which comes out from the extrapolation of our local 
(solar) coordinate system, always exceeds the actual physical volume (due to 
the presence of wormholes). The ratio Veoora/Vpnys = Q defines the porosity 
coefficient of the foamed space. When we use the extrapolated coordinates we 
always overestimate (by the use of the Gauss divergence theorem) the actual 
intensity of a source of gravity or of an incoherent radiation. In gravity, the effect 
displays itself as the presence of Dark Matter. Hence, the porosity coefficient of 
the foamed space Q can be related to the ratio of Dark Matter density to the 
density of baryons in the Universe, i.e. Q = Qpy/Qy. Analogously, the same 
relation holds true for the ratio of two components of radiation (diffuse phone 
and discrete sources), i.e. Q = OQaiffuse/Qdiscrete. The relation 


Qpm/Qb x Qdif fuse /Qdiscrete 


is the basic indication of a geometrical (topological) nature of DM effects. 


2 Random Topology of Space 


In order to set a general frame for the study of a foamed space, let us start with 
a toy example where the space is a cylinder of radius R. The metric is the same 
as for the standard flat Friedman model 


ds? = dt? — a? (t)(dz? + dy? + dz’), (6) 


but one of the coordinates, say z, is periodic (z +27 R = z). In what follows, for 
the sake of simplicity we neglect the dependence of the scale factor on time in 


(6), i.e. consider the Minkowsky space as the coordinate space. Thus the actual 
values of the coordinate z run through the fundamental region z € [0,27R]. 
Such space can be equally viewed as a portion of the ordinary R? between two 
plane mirrors (at the positions z = 0 and z = 27R). An observer, who lives 
in such space, may use the extrapolated reference system (i.e., z E (—o0, 00)), 
however he/she easily notices that all physical fields are periodic in z. Consider 
the Newton’s potential ¢ for a point mass M. In this space the exact expression 
can be easily found from the standard Newton’s potential by means of the image 
method. Indeed, the periodicity in z means that instead of a single point mass 
Mo6(r —r') at the point r’ the observer will actually see an infinite series of 
images 


i(r—r!) > K (rr!) = Y -2 + 20Rn)6(w@——2')5(y-y'), (7) 


n=—Cco 


and the Newton’s potential for a point source at r’ = 0 takes the form 


p=-GM X` 1/\/p? + (z+ 20Rn)’. (8) 


n=— oo 


On scales r < R we may retain only one term with n = 0 and obtain the 
standard Newton’s potential for a point mass ¢ ~ —GM/r, while for larger 
scales r >> R the compactification of one dimension will result in the crossover 
of the potential to ọ ~ cM inr (note that this is indeed the shape of the 
potential that one reads from the observed galaxy rotation curves). 

The anomalous behavior of gravity indicates that DM effects show up at this 
model on the scale of distances of order R. Indeed, let us consider a box of the 
size L and evaluate the total dynamical mass within the box 


L 
Mot (L) = M „E (r,0)dV = M (1 + isa) . (9) 


Thus, if the observer is ignorant about the real topological structure of space he 
should conclude the presence of some extra (odd) matter. The reason is obvious, 
when we expand the coordinate volume it covers the physical (or fundamental) 
region of space many times and we respectively many times account for the 
same source (i.e., images of the actual source). Thus the increase of the total 
mass is fictitious. In the simplistic model under consideration the volume of the 
fundamental (physical) region behaves as Vphys = L? for L < R and Vphys = 
2rRL? for L > R. We note that at large distances L > R the parameter 
Q (L) = Mtot(L)/M —1 can be used to estimate the actual value of the physical 
volume: Vpnys(L) = L3/Q (L), i.e., Q is the “porosity coefficient” of space at 
scales L ~ R. 

The space discussed above is rather simple: for an extended source we will 
see a countable set of its images without distortion. Therefore, one can easily 
detect the fundamental region of space and avoid consideration of fictitious 


sources. In the case of a general foamed topological structure this is hardly 
possible. Nevertheless, whatever the topological structure of the manifold is, we 
can apply the method of images: every topology can be achieved by introducing 
a certain equivalence relation in R? and gluing equivalent points together. Thus, 
a space of non-trivial topology is completely defined by indicating for every point 
r! € R? the set E(r’) = {fi(r’), fo(r’),...} of the points equivalent to it. In other 
words, a point source at a point r’ in the fundamental region is accompanied 
by a countable set of images, or “ghost” sources: 


i(r-r')4+K(nr)=d(r-r')+ YO ar- flr’) (10) 


fir )EE(r’) 


where f; (r’) is the position of the i-th image of the source. 
For example, consider any source for radiation J(r,t). Then according to 
(10) the electromagnetic potential A (0, t) is described by the retarded potentials 


J Ir|/e J r)|/e 
A= - [|= = E dV + a oA TA (11) 


The first term of this formula corresponds to the standard, “direct” signal from 
the source, while the sum describes the multiple scattering on the topological 
structure of space. A similar formula is obtained for the gravitational field. 

It is clear that all physical Green functions for all particles acquire the same 


structure 
Grora(0,7) = Go (0,r)+ XO Gi(0,filr))- (12) 
fi(r)EE(r) 
Formally, one can use the standard Green functions, while the scattering will 
be described by the bias of sources 


Jtotat (r,t) = J (r,t) + J b (r,r) J (r’, t) dr, (13) 


where b (r,r) = K (r,r')—ô(r— r’), i.e. we excluded the actual point source. In 
gravity the second term in corresponds to the DM contribution (e.g., see 
[4]). We note that in general the bias b(r,r’) is an arbitrary function of both 
arguments, which means that the nontrivial topological structure is capable of 
fitting an arbitrary distribution of Dark Matter. 

The function K(r,r’) unambiguously defines the topological structure of the 
physical space. However, for a general foamed structure of space (a gas of 
wormholes) this function has a quite irregular character, i.e. it is not directly 
observable. One has to introduce a measure on the space of all 3D-manifolds 
and average the function K over this measure. The resulting function 


K(r,r’) = 6(r — r’) +. b(r, 7’) 


gives the (average) density, at the point r, of the images of the point r’. 


Because of the averaging, the irregularities are smoothed out, hence the bias 
function b(r,r’) is observable. Indeed, the averaging of and gives 


Gtotal (0, 7’) = G (0, r’) + [07 r')G (0,r)) dr (14) 
for Green functions, and 


Ptotal (r,t) = p(r,t) + I b (r, r") p (ri, t) dr’ (15) 


for the density of matter. Therefore, when we can distinguish two components in 
the observed picture of the distribution of, say, gravity sources: discrete sources 
and a diffuse phone, the discrete sources can be identified with the first term 
in the right-hand side of (Q5), i.e. with “actually existing” sources, while the 
diffuse halo can be identifies with the second term, “the images”. Then, by 
comparing the observed distribution p(r’) of actual (discrete) sources with the 
observed DM distribution 


matie J B (r, r’) o(r)dr', (16) 


one can extract an information about the structure of the bias b. In fact, the 
homogeneity of the Universe requires from b to be a function of (r — r’) only 
(which means that the form of DM halos does not, in general, depend on the 
position in space). In this case, the Fourier transform of (16) gives 


Phaio(k) = b(k)p(k), (17) 


which defines b uniquely. As we show in the next Section, the bias b extracted 
from the DM observations in this way has both a very simple form and a trans- 
parent theoretical meaning. 

Note that being an averaged characteristics, the bias b does not determine 
the topology of space completely. Along with the one-point distribution K (r, r’), 
one can consider joint distributions of images for several sources: 


Krlri, see Tn; Tris PAR Th) 
which is the averaged density of the images of the points r},...,7/, at the points 
T1,--++,;T- Only when all the functions Kn, n = 1,2,..., are determined, one will 


have a full description of the structure of the foamed physical space. However, 
the one-point bias functions 6(r,r’) carries the most important information. 

Thus, consider a source of radiation, constantly emitting light with the fre- 
quency w, i.e. we have a density of the EM current J(r’)e™* such that 


(I(r) I* (r5)) Z ôlri = 15) Tsource(T1), (18) 


where Isource(7) is the spatial distribution of the intensity of the source. In 
order to take into account the effects of the non-trivial topology of space, J(r) 
should be modified according to , Le. J(r1)J*(r2) transforms into 


fx (r1, ri) (rate) J(r1)J* (r3) Jdr! r= | K( ri, f K (r2, r) Isource(r dr. 


Averaging over different topologies gives 
(J(r1)J*(r2))iotal = fk (r1, T2; r, r’) Line ar; (19) 


where K2(r1, r2; r’,7’) is, by definition, the joint distribution of a pair of images 
of the point r’. 

The points rı and rə can be images of the same point r’ if and only if 
they are images of each other. Therefore, K2(r1,1r2; 1’,r’) is proportional to 
K(r1,72) = 6(r1 — r2) + b(r1, r2); more precisely 


Kə (r1, r2; fs r’) = òlr = r2ə)K (r1, r") a b(ra, r2)P(r1,72,1') (20) 


where we denote as P(r1,1r2,r’) the density at the point r2 of the distribution 
of images of the point r’ under the condition that the point rı Æ r2 is an image 
of T2. 

As we see from (19), (20), while the phases of the source current J(r’) are 
delta-correlated (see (18)), there appear long-range correlations in the density 
of the total current — due to the term proportional to b(r1,r2) in the kernel 
Ko. However, the characteristic wave length in b(rı — r2) is of order of galaxy 
size, i.e. it is unimaginably larger than the wave length cw of the light emitted. 
Therefore, the contribution of the coherent part of the total current to the 
radiation is completely negligible: by (19),(20) we find 


(J(r1)J*(r2))sotar = O(71 — r2) / K (r1,7’) Isource(r’)d?r’ + long wave terms, 


which gives the following formula for the total intensity of sources (actual plus 
ghost ones) 


Ttotat (r) = JE (r, r’) Isourcelr er = TL eunéelh) + / b (r, r’) L source (r') dr’. 
(21) 
Comparing with (IG), we see that the distribution of a diffuse radiation phone 
associated to a luminous source coincides with the distribution of dark matter 
in the halo of the same source. 

Note that for a non-stationary remote source of radiation the picture is more 
complicated. A momentary pulse at some point will create a spherical EM wave 
emanating from the point — and from its images. On the front of the wave only 
a small number of images will give an essential contribution, namely those which 
have comparable and shortest optical paths. This will lead to an interference 
picture on the front. We note that due to wormholes the signal from some 
images can reach an observer even earlier than the basic signal. Only with time 
elapsed, as the larger and larger number of images contribute, the interference 
picture disappears, and the diffuse radiation phone given by (21) establishes. 

In conclusion of this section, we recall that the observed homogeneity and 
isotropy of space require from the topological bias b(r, r’) that defines both the 


DM distribution (16) and the distribution (21) of the sources of diffuse radiation 
to be the function of the distance |r—r’| only: b(r,r’) = b(|r—r’|). The integral 


Q(L) = 4n [ R?b(R) dR (22) 


characterizes then the distortion of the coordinate volume or the porosity of 
space (i.e., 1/Q gives the portion of the fundamental region or the volume of 
the actual physical space in a coordinate ball of the radius L). In general there 
can be both a situation where Q(L) tends to a finite limit as L > oo and then 
Q(ov) defines the total amount of DM (Q = Qpm/Qb = Qiiffuse/ discrete) 
and the case where Q is unbounded. The last case indicates the presence of a 
certain dimension reduction of space at large distances (e.g. when Q(L) ~ L® 
the dimension of the physical space reduces to D = 3 — a). 


3 Topological bias: empirical and theoretical ap- 
proach 


In this Section we derive a formula for the bias function b(|r—r’|) and show that 
it fits the observed picture of DM distribution quite well. While in empirical 
considerations it is more convenient to view b(R) as a bias of sources (which 
means exploring the laws (6) and (21)), we achieve more theoretical insight 
when choose an equivalent description of the random topological structure of 
space by means of the bias of Green functions (see (14)). This means that instead 
of saying that each material point is accompanied by an infinite set of images, we 
say that each source excites an infinity of fields. Indeed, on a connected manifold 
of non-trivial topology there is an infinite number of geodesics connecting any 
two points. So the light emitted at a point P arrives at a point Q by an infinite 
number of non-homotopic ways. We may associate a separate EM field with 
each homotopy class: each of the fields propagates independently, but they sum 
up when interact with matter. When we describe things in R? by means of 
the bias functions, we thus associate a separate field to each term in the right- 
hand side of (12). These terms differ by positions of the images f;(r). In our 
picture, where the topology is random, there is no preferred position for the 
i-th image, hence we have a system of an infinite number of fields {A;} which is 
symmetric with respect to any permutation of them (in other words, the fields 
are identical). 

It is widely believed that the effects of quantum gravity should lead to a 
cut-off at large wave numbers. The cut-off at A means that the photons with 
wave numbers |k| > A are never excited. We say that the field does not exist 
at such k. One can describe a cut-off of a more general form, by introducing 
a characteristic function x(k): at x(k) = 1 the field with the wave number k 
exists, while at x(k) = 0 it does not. Because of the renormalizability of all 
physical field theories, the question of the determining exact form of the cut-off 
of a given field is of little importance. However, for the system of an infinite 
number of identical fields {A;} the cut-off function acquires a meaning. 


Indeed, let us define N(k) = X; x:(k) where the sum is taken over all the 
fields A;. Thus, N(k) is the number density of fields which exist (i.e. which 
are not forbidden to create particles) at the given wave number k. Here, the 
existence of the cut-off means that N(k) can be finite for all k. As the fields 
sum up when interacting with the matter, the values of N(k) greater than 1 
lead to a stronger interaction than in the case of a single field. For example, 
consider a Newtonian potential 


Ad = 4ryp. 
In the Fourier representation we have 


— Ary 


o(k) = p(k), (23) 


If there exist N(k) identical Newtonian gravity fields with the wave number 
k, each of them satisfies (23), while the effective potential (that which acts on 


matter) is given by dep f(k) = ye ) i(k) and satisfies, therefore, 


desp(h) = SE N(R) ph). 


This is equivalent to a renormalization of the source density 
p(k) > N(k)p(k), 

and comparing with gives 
N(k) — 1 =b(k). 


Thus, the Fourier transform b(k) of the topological bias function can be inter- 
preted as the excessive number density of fields (gravity or EM) at the wave 
number k, i.e. it is determined via a cut-off function. 

Although the problem of determining the exact shape of the cut-off is usually 
considered hopeless because the full quantum gravity theory has not been devel- 
oped, an approach developed in [2] allows one to derive possible types of cut-off 
by means of simple thermodynamical models. For example, assume that the en- 
ergy density and the total excessive number density of fields N = [(N(k)—1)d?k 
are finite. We also assume that M is a conserved quantity (along with the en- 
ergy). Then the shape of the function N(k) is determined uniquely by the 
condition that the system of the identical free fields is in the thermodynami- 
cal equilibrium (one should only choose the statistics for the fields and fix the 
values of thermodynamical parameters). Indeed, the state of the system with 
N (k) identical free fields at the wave number k is determined by the numbers 
ni(k), i = 1,..., N(k) of the particles with the wave number k for each field. In 
the case of Fermi statistics for the fields (that has nothing to do with the statis- 
tics for the particles which remains Bose), there cannot be more than one field in 
the given state, i.e. for every given k all the numbers n;(k) should be different. 
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The energy density at the wave number k equals to wk pally ) ni(k), where wk 
is the energy of a single particle; as we deal here with massless fields, we take 
wp = |k| (we put h = c= 1). In what follows we assume Fermi statistics for the 
fields (Bose statistics leads to a similar result [2] [5], however the computations 
in Fermi case are simpler). Then, the state of the lowest possible energy (“the 
ground state”) corresponds to {n1(k),... nn (k)} = {0,1,...,N(k) — 1}. 
This gives us the energy |k|N(k)(N(k) — 1)/2 at the wave number k. The 
total energy density is thus given by f ILN (k) (N (k) — 1)d?k. The ground state 
corresponds to the minimum of the total energy density. As the total excessive 
number density of fields M = f (N(k) — 1)d?k is assumed to be conserved, the 
problem of finding N (k) reduces to minimizing f |k|N(k)(N(k) — 1)d?k under 
the constraint f (N(k)— 1)d?k = constant. This gives us 


N(k)=1+ hl , 


where the “chemical potential” u is fixed by the value of M. For the bias function 
b this gives 


H 

7 — for |k| <p, 

w= | ao N (24) 
0 for |k| > yp. 


One can make different assumptions and, perhaps, arrive at different for- 
mulas for the bias. However, this simplest bias function provides a very good 
description of the observed distribution of DM. Indeed, in the coordinate rep- 
resentation bias takes the form 


Bi) = af Own) BOS (25) 
0 
= a (1 — cos (ur)). 


As it was shown in [6], by choosing u = 7/(2Ro) where Ro is of order of a 
galaxy size (i.e. a few Kpc), bias applied to spiral galaxies produces the 
pseudo-isothermal DM halo p = po R2,/(R2 + r°), where Rc is the core radius 
which has the order of the optical disk radius Ro ~ Ropt. We note that this is 
in a very good agreement with the observations (see [8]). In fact, by fitting one 
parameter p in accordance to Tully-Fisher law [14], relations and quite 
accurately represent the whole variety of the observed galaxy rotation curves 
[6] (we recall that bias is derived from thermodynamical considerations, 
so it is quite natural to allow the chemical potential u fluctuate in space; exact 
mechanisms governing these fluctuations are described in [4] [6]). 

From (22),(25) one can find that starting with the galaxy scale the porosity 
of space behaves as Q (r) ~ r/Ro. Thus the total dynamical mass for a point 
source within the radius r increases also as 


M (r) ~ M (1 +r/Ro). (26) 
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Importantly (see the previous Section), the same conclusion holds for the lu- 
minosity of the point source (i.e., for a galaxy or an X-ray source). There- 
fore, one can not immediately conclude from (26) a linear growth of the ratio 
Mtot(r)/My(r) of gravitational (dynamical or lensing) to the barionic mass: the 
result depends on how much of diffuse radiation is discarded at the observations. 

Observations suggest that the number of baryons within the radius r behaves 
as M (r) ~ rP? with D ~ 2 (see e.g. Refs. where the ~ r? behavior 
was reported up to at least 200 Mpc). Thus, the observed barionic density Qs 
falls inverse proportionally to the deepness of the observations and is well below 
1. In the standard picture the total gravitational mass grows as ~ R?, as it 
should be in a homogeneous Universe, so the linear growth of Mtot(r)/Mi(r) 
predicted by bias is indeed consistent with observations. However, the 
linear growth starts to show up with the scales larger than cluster size, while 
the reported mass to luminosity ratio remains approximately the same on the 
galaxy scale and on the cluster scale. To resolve the problem, we invoke the 
results of [I] where there was demonstrated that the intracluster gas clouds may 
not carry dark matter. In our picture this is indeed the case, as the intracluster 
cloud is an extended source of X-ray radiation, of size much larger than Ro. 
Thus, the associated diffuse phone sums up with the “direct” signal, so all the 
ghost sources of gravity that lie within the cloud are visible as well. This means 
the absence of “dark” matter in the cloud or, in other words, that the number 
of barions in the cloud is greatly overestimated — most of the contribution to 
the cloud luminosity is given by the diffuse halo, i.e. by fictitious sources due to 
the non-trivial topology of space. It is easy to check that correcting the barion 
density of the intracluster gas in accordance with (21), provides indeed the 
linear growth of Mtot(r)/Mo(r) starting right from the galaxy scale. 

Note that at very large scales the diffuse radiation can hardly be separated 
from the very faint sources. Therefore, the picture of the homogeneous distribu- 
tion of matter (i.e., of the Friedman Universe) is restored. In fact, an arbitrary 
foam-like structure of space (i.e., any choice of the bias b(r)) agrees perfectly 
with the observational large-scale homogeneity and isotropy of the Friedman 
Universe provided that the actual physical volume Vpnys (r) = 4/37r3/Q (r) 
(the volume of the fundamental region of the coordinate space) is homogeneously 
filled with matter. Indeed, in this case the number of actual sources within the 
radius r behaves as the physical volume Np (r) ~ Vpnys (r) ~ r?/Q (r). Along 
with the actual sources we always observe images (DM and diffuse radiation) 
and every source produces AN ~ Q (r) additional images. Thus the total num- 
ber of images behaves always as M (r) Q (r) ~ r3, i.e., produces a homogeneous 
distribution. 


4 Conclusion 
In conclusion, we briefly repeat basic results. First of all the concept of space- 


time foam introduced by Wheeler can be crucial in explaining properties of the 
present day Universe. The random (“foamed”) topological structure leads to the 
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fact that every discrete source in the sky should be surrounded with a specific 
halo (a random distribution of images). We call this phenomenon a topological 
bias of sources. In gravity such halo modifies the standard Newton’s law and 
appears as the Dark Matter phenomenon. In particular, the Universal rotation 
curve (URC) constructed in [6] on the basis of the topological bias shows a very 
good fit to the empirical URC [8]. We stress that in a general foamed space 
the bias b(r,r’) is a random function of both arguments which means that the 
form of the DM halo can arbitrary vary in space. By other words any observed 
distribution of DM can be easily fitted by a proper choice of the foamed struc- 
ture. However, the simplest bias function which we derived theoretically from a 
basic physical (thermodynamical) considerations seems to give a quite accurate 
account of the DM effects in a huge range of spatial scales. 

As it was demonstrated in this paper, in the foamed space the halos around 
discrete sources are actually not dark, but form the diffuse background of ra- 
diation. Moreover, the ratio of the two components (the diffuse background 
and discrete sources) is exactly the same as the ratio of DM and baryons 
(Qpu/Qv = Qdif fuse/Qdiscrete): 

We note that the foamed picture of our Universe allows to explain the prob- 
lem of missing baryons. Recall that the direct count of the number of baryons 
gives a very small value Q, ~ 0.003 for the whole nearby Universe out to the 
radius ~ 300hs5 M pe e.g., see [9]. In our picture, this means only that at the 
radius ~ 300h55 M pc the actual volume is ten times smaller, than in the Fried- 
man space (Vpnys ~ 0.1Vr), i.e. the actual density is ten times bigger which 
reconciles the observed small barion density with the primordial nucleosynthesis 
constraints. 

We stress that any homogeneously filled with matter foamed space (i.e., an 
arbitrary choice of the bias function b(r,r’)) agrees perfectly with homogeneity 
and isotropy of the Universe and does not contradict to the standard Friedman 
model. The general foamed Universe can be viewed as the standard Friedman 
space filled with a gas wormholes. In such a picture the Large Scale Structure 
has an equilibrium character, for it reflects the foamed topological structure of 
space (i.e., the distribution of wormholes) formed during the quantum period of 
the evolution of the Universe. 

Finally, we have demonstrated that in a foamed space any non-stationary 
and sufficiently remote signal is accompanied with a formation of a specific in- 
terference picture at the front of the wave (stochastic interference) which rapidly 
decays. 
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